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Abstract 

We present a steepest descent calculation of the Kallen-Lehmann spectral density of two-point functions involving 
complex conjugate masses in Euclidean space. This problem occurs in studies of (gauge) theories with Gribov-like 
propagators. As the presence of complex masses and the use of Euclidean space brings the theory outside of the strict 
validity of the Cutkosky cut rules, we discuss an alternative method based on the Widder inversion operator of the 
Stieltjes transformation. It turns out that the results coincide with those obtained by naively applying the cut rules. We 
also point out the potential usefulness of the Stieltjes (inversion) formalism when non-standard propagators are used, in 
which case cut rules are not available at all. 



1 Introduction 

In the setting of local quantum field theory, all information on the spectrum resides in the Green functions of the 
theory, as they build up the scattering 5-matrix. A special class of Green functions are the two-point functions. 
We should not only consider the elementary propagators of the theory, but also the two-point functions of local 
composite operators, which correspond to propagating bound states, are relevant. 

Scattering amplitudes are in general analytic functions of the Lorentz invariants constructed from the external 
momenta characterizing a given process. In the present work we shall only be interested in two-point functions 
which can be considered as functions of the external momentum k"^ = s, analytically continued to the complex 
s-plane with a cut on the real axis starting at s = tq > 0, which is identified as the threshold for multi-particle 
production. Such a two-point correlation function can be cast in the form of a dispersion relation, that is, an integral 
representation written in terms of the function's discontinuities in the complex plane. For the two-point function 
of an arbitrary operator O we can write the general expression 

{o{k)oi-k)) ^ j-(s) = - r dT ^^"lzm . (1.1) 

The imaginary part in the argument of the integral stands for the discontinuity of the function J-"(s) across the 
cut. This ought to be a positive definite quantity proportional to the total cross-section, as demanded by the optical 
theorem. The representation is known as the Kallen-Lehmann spectral representation. For a concise treatment, 
let us refer to e.g. HI E) [3] HJ . 

This Kallen-Lehmann spectral representation can be constructed, for the case of real masses in Minkowski space, 
through the use of the Cutkosky cut rules fT '5], which state that, in order to calculate the discontinuity associated 
with a Feynman diagram, we just need to cut propagator lines in the diagram in all possible ways, thereby replacing 
the corresponding propagators by (5-functions localizing to the physical phase space, and then to sum over all these 
contributions. Physically, we place the particles running in these lines on-shell. 
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Specifying for a moment to the strong interaction, it is well-known that QCD exhibits confinement, i.e. the color- 
charged elementary degrees of freedom (quarks and gluons) do not propagate in the physical world. The observ- 
ables are the color-neutral bound states, corresponding to the gauge invariant operators constructed from the gluon 
and quarks fields, being the mesons, baryons, glueballs and quark-gluon hybrids. A big contemporary challenge 
is to extract information on that nonperturbative spectrum using theoretical, phenomenological and/or numerical 
lattice tools. In the so-called sum rules approach to the QCD spectrum, spectral representations play a pivotal role, 
we refer to the comprehensive tome [6], to [TuSJ for the seminal works, or to |4| for a pedestrian's introduction. 

In many approaches to QCD, including the sum rules, an important role is played by the elementary propagators, in 
particular the gluon propagator, which is usually studied in the Landau gauge, d^A^^ = 0. The elementary propaga- 
tors are the fundamental building blocks of more complicated Green functions, including the ones corresponding 
to bound states 191 [TOl [TTI . For some recent advancements on the gluon propagator, including numerical lattice 
results, let us refer to the non-extensive list of lfT2l[T3l[T4l[T5l[T6l[T7l[T8iri9, 20 2l1l22ll23ll24ll25ll26ll27ll28l. 



It is apparent that the Landau gauge gluon propagator exhibits a violation of positivity 12911301 . a fact which also 
received direct lattice confirmation Il3ni32l . This violation of positivity means nothing more than that there cannot 
exist a Kallen-Lehmann representation fTT^, hence the Landau gauge gluon is unphysical, and it can be interpreted 
as a reflection of confinement. In 1.18. .19J . it was discussed how 

can qualitatively and quantitatively describe the lattice predictions for the (infrared) Landau gauge gluon propa- 
gator, and this up to fc ~ 1.5 GeV. The expression ( |L2| l can be seen as a dynamically improved version of the 
so-called Gribov propagator, 

^(fc^)-^^ (1-3) 

which attracted a lot of attention as a prototype of a confining gluon propagator, and which is related to the dealing 
with the issue of gauge (Gribov) copies in the Landau gauge Il33l[34ir35i . Given that ( |L2| l works out well, it might 
subsequently be used in explicit computations, for example when one would like to learn something on the glueball 
or meson spectrum. 

We notice that the propagator jl.l) can be decomposed into 2 Yukawa propagators, with either 2 real, or 2 complex 
conjugate masses, depending on the relative size of the scales present. The case of 2 complex conjugate masses 
seems to be realized [19], and a (tree level) positivity violation is then guaranteed. The unphysical degrees of 
freedom corresponding to those complex conjugate masses were called i-particles [36|. Also the so-called Stingl 
propagator ll37l[38l . 

Dik')= , f,, , (1.4) 



fc4 + ,n2/s2 + X4 

can exhibit a pair of complex conjugate poles, as a special case of ( |1.2| i. 

If we are interested in a spectral function inspired approach to e.g. bound states in theories with propagators of the 
type ( |1.2| l, we of course would need to know how to compute with such a propagator. Research in this direction 
was initiated with the work on a toy model [361, where the necessary spectral densities were computed by clever 
manipulations, which unfortunately do not have general applicability, or, more realistically, for gauge theories in 
||39l l40| by using the Cutkosky cut rules and assuming that the results can be continued to Euclidean space and the 
case of complex conjugate masses. 

The main purpose of this paper is to show that the cut rules results, which we summarized in Section 2, are 
actually correct. We shall therefore bring to attention in Section 3 the Stieltjes integral transformation, as well as 
its inverse, as studied by Widder |41 1, which allows to compute directly in Euclidean space the spectral density, 
even in cases when there is no cut rules prescription. We shall work out an introductory example in Section 4. In 
Section 5, we comment on the use of subtracted spectral representations to avoid divergences. We then present 
in Section 6 in full detail the most interesting case of 2 z-particle propagators with complex conjugate masses, 
wherefore a completely analytical treatment can also be given. In our concluding remarks. Section 7, we point 
out that the Widder result to compute the spectral function could perhaps also be used in a numerical context 
when more complicated propagators than e.g. ( |1.2| i are used, which are completely inadequate for fully analytical 
computations. 
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2 Using the Cutkosky cut rules 



2.1 Preliminaries 

We assume that O is a generic (composite) operator containing 2 fields. We are then interested in the study of the 
following 1-loop expression 



where /(p, k—p) is a function of the Lorentz invariants constructed from the momenta flowing in the loop. For sim- 
plicity, we may assume O to be a scalar function, thus a fortiori also /(p, k — p). If O would carry Lorentz indices, 
say O = 0^i...p„, the eventual correlation function (O^i...^,, = J^pi...^„,.i...,.„ (fc, toi, mz) 

would be a tensor object. One can always decompose this tensor into a suitably constructed basis (using the vector 
k^ and trivial unit tensors), and by projecting on the basic tensor structures, one can construct the spectral densities 
associated to each structure. The latter problem can then be reduced to that of the scalar case by contractions. For 
example, assume that the two-point function is a rank 2 tensor, J^^^[k). Each rank two tensor can be decomposed 
into its transversal and longitudinal part, according to 

^^.(fc) = ^T(fc) - ^) +^L(fc) , (2.2) 



SO that knowledge of Fxik), resp. J^L{k), follows from the (scalar) quantity P^u{k)J-^u{k), resp. L^^{k)F^,y{k). 
Using these scalar objects, one can the construct the spectral densities corresponding to the transversal and lon- 
gitudinal piece of the two-point functionb. The here described procedure can be generaUzed to more complicated 
tensorial structures. 

We thus wish to find the spectral representation of the integral \2.\\ , and we shall therefore use the Cutkosky 
cut rules. In order to do so, we must assume that m\ and are real, and that we are in Minkowski space time. 

By cutting the lines comprising this 1-loop diagram and, according to the rules, performing the replacement Em 

1 



(p2 _ rnl) 



27r0(p°)5(p2 _ ml) , (2.3) 



the spectral function is proportional to the discontinuity, p(t) oc Disc T{t) — 2Im F[t), and it can be determined 
from 

lmF{k,m,,m2) = \ J ^ [(27r)^^^((fc - p)")<5((A: - p)^ - m?)0(pO)<5(/ - m^)/^ fc - p)] . (2.4) 



In order to evaluate (2.4i we will work in a frame where k^ = (fc°, 0) = {E, 0). We have then 



Im :FiE') = \j j^^z^ [OiE - p")5{iE - p^ - ^l,)e{p^)5{{py - ujl,)f{p, k-p)], (2.5) 

where Wp.^ = yj-j? + mf. In the evaluation of this integral one must observe that the function /(p, k—p) comes 
from possible derivatives in the operator 0{k) and that it is a Lorentz scalar It can thus only be a function of E, 
nil and 1712, which comes from the scalars p^ = mj, (fc — p)^ = nif and p ■ {k — p) — \{E^ — m\ — m^). After 
the evaluation we can write the Lorentz invariant spectral function p{t) by going to an arbitrary frame, whilst 
replacing E"^ — > t. 

We are ultimately interested, however, in the case where the masses involved are complex and the momenta are in 
Euclidean space. We then consider expressions like \2.5\ as functions of real masses and Minkowski momenta and 
will analytically continu^it to complex mass values in Euclidean momentum space. We want to draw attention 
here that, if we would immediately assume ra\ and rri^ to be complex, we cannot apply the cut rules. We cannot 
put an z-particle physically on-shell, or mathematically spoken, ( |2.3| l does not make much sense for m? G C while 
p2 g M. 



1 . A certain care is needed when speaking about "analytic continuation", as we are considering a function depending on two variables, viz. 



and m| . 
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2.2 Naive application of the cut rules 

We will now illustrate the procedure just described by computing spectral functions in different dimensionalities 
for the two-point function ( |2.1| i with f{p,k — p) — 1. This corresponds to the toy model introduced in ll36ll . We 
are thus interested in finding a spectral representation of the form 

(2.6) 

for the following Euclidean two-point function 

/d'^p 1 1 
(2^0^ p-^+mj {p - kY + ml ' ^^-^^ 

with conjugate masses parametrized as 

rnf — a + ib , m\ — a ^ ib . (2.8) 
We assume thaj^a > 0, and without loss of generality we may take & > 0. 



We shall therefore first use the Minkowski formula ( |2.5| ), and integrate over in ( 2.5 1, obtaining 

1 /■ S'^^^^v 1 

HE') - 2 j ^ [^(^ - --^)^((^ - -^-f - <^)\ • (2.9) 

The integrand depends only on so that we can write 



d-i 



where we also took into account the theta function. The remaining delta function can be cast in a more convenient 
form using the property (5 = |7^^(P^-Po) = 2|po||g'(p!)l ^^^1 ~ where is such that g(pg) = 0. 

In our case, p^ satisfies 



g{f^) = E-^f^ + ml-^pl+ml = Q. (2.11) 
The integral ( |2.10| i can now be readily evaluated 

^ ' (27r)('i-2) r(^^) 4£; ' ^ ^ 



bol = ^ ^ , (2.13) 



where 



Having done this, we can now pass to the Euclidean case with the masses given by ( |2.8| l, by evaluating this 
expression ( |2.13| l, 

Ipo I = 7^ . (2.14) 

The expression ( |2.12| i can then worked out further for the various dimensions. For (i = 2 we obtain 

Im Fd-2{E'') = ^ ■ (2.15) 



We eventually conclude that 



Pd=2{T) = —-^====== , (2.16) 

2tt - 452 - 4ar 



2. This in order to avoid tachyon instabilities. For 6 = 0, this is immediately clear. In the context of the model of 1 181 . a itself also corresponds 
to the mass of another degree of freedom, as such a is supposed to be positive. 
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using E'^ T and the equivalence pij) = ilmJ^(r). 
Analogously, for d = 3 and d = 4, we find 



lm^..3(i^^)^i^^p..3(r) = ^-^, (2.17) 



The threshold tq is in all cases given by 

To = (mi + 7712^ = 2 (a + ^/a^ + b^^ . (2.19) 



3 Survey of the Widder Stieltjes inversion operator 



In this section, we shall first refresh the concept of the Stieltjes integral transformation, and then discuss how the 
inverse transformation can be found. 



Let us assume a function F{x), defined for a; > 0, and implicitly define another function p{t) by means of 



Fix) 



t + X 



(3.1) 



where we assume that the latter integral exists. This operation defines the Stieltjes integral transformation 11411 l42l . 
Upon comparing and it is clear that a Stieltjes representation of the function F{x) is nothing else 
than the Kallen-Lehmann representation if F{x) would be a two-point function in Euclidean space. If p{t) — for 
t < A, A sets the threshold tq. 



The spectral density p{t) can be reconstructed from Il4ni42ll 

p{t)= lim (-1)"+^^ 



(3.2) 



Let us first prove this statement along the lines of 141], as the tools of the proof shall turn out to be useful for our 
later analysis too. We start with 



t + U 



n + 1 







2n+l 



n + 1 



-j-oo 



dup{u)u"-dl- 



(i + u)«+2 



(3.3) 



It seems to be more involved to compute the latter derivative, but we can use a nice trick to do so BTI . Defining 



g(t,u) 



(3.4) 



+ m)"+2 ' 

we notice that this is a homogenous function of order —1, as for any ^ > 0, 

g{itju)^e-'g{t,u). (3.5) 
Deriving w.r.t. £ and setting £ — 1 leads to the Euler characterization of homogenous functions. Specifically, 

tdtg + udug = -g, (3.6) 

or for t > 0, 



dtg -du i^g^ 



(3.7) 
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If g(t, u) is homogenous of order —1, so is f u), hence we can iteratively employ p.7[ ) to compute 

^/ -I'.K.K+lort 1 ^ ,rt+l (2"- + 1 /TON 

^ ^ "(i + w)"+2 (n + 1)! (t + M)2«+2 ■ ^■'■"^ 



We are thus lead to 



u 



with 

Mu,t)=ln^— (3.10) 
In the limit n — > cx), the latter integral is ideally suited for a steepest descent evaluation, since 

^nM",i)L=t = 0, dlh{u,t)\^^^^-^<Q. (3.11) 

Doing so, we find 

n->+oo (?l!)2 _/q u n-i.+oo (n!)^ <; _/g 

(3.12) 

with 

h{t,t) ^ -2\ii2 , (3.13) 

and, for n +oo, 

due-("+i)^T^ ^ / due"^"+^)^i^ ^2 / ^. (3.14) 
J-oo V^ + n 

Recalling Stirling's formula 

n! -> V27rn (^-j , for n +00 , (3.15) 

we can write that 

n2^^' ^ -j^4"+iy^, for n-^ +00. (3.16) 



We then find 



1^+00 (n!)2 t Jo 



hereby proving p.2| l. 

Once having determined p{t) via p.2| i, one can then define 



(3.18) 



for z e C, which is analytic, with the exception of a branch cut on the negative real axis, z < [42]. Obviously, 
there will be no branch cut for z G [—62, —Si] if p{t) = for i e [61, 62], with 61 > 62 > 0. Using Cauchy's 
formula, one can then see that it actually holds that 

p(t) ^ lim \F(-t - ie) - F(-t + ie)] , (3.19) 

which corresponds to the discontinuity of F{z) along the negative real axis. En route, this explains the dispersion 
relation ( |l.l| l. 

Apparently, p.l9| l would provide us a with a much easier way to compute the spectral density p{t), but the difficulty 
is that in most cases, we do not know how to evaluate F{z) from its original (integral) definition for z ^ Ej. In 
contrast with this, p.2| i only needs F{z) for z e ■ This is exactly what we need, since we can always evaluate 
the Euclidean momentum integrals defining a two-point function F{k'^), in which case we can assume G Rq . 
Afterwards, F{z) is defined by means of p.l8| l for all z E C. 

In principle, p.l8| l also makes sense when p{t) is a distribution. A basic example is p{t) — d{t — m?), which leads 
to the familiar Yukawa-propagator F{k^) = ^,2 . which evidently has no branch cut. 
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4 A first application: two (positive) real masses 

As a warming up exercise, we shall treat here a well-known textbook example, and compute the spectral density in 
the case of 2 real masses. This was already treated in [3_|, albeit in Minkowski space. We shall work in Euclidean 
space, and we start from (|2.7|i. We introduce a Feynman parameter x, yielding 



(4.2) 



The substitution q = p — xk gives 
with 

= xk^ + xm\ +ml- xmj - x^k^ . (4.3) 
Usage of the standard Euclidean space formula 

A 1 _ 1 ,.2^d/2-„r(n-d/2) 



(A^)"/^-" ' 7 ' (4.4) 



(27r)rf (g2 + A2)« {AttY/^^ ' r(n) 
leads to 

^'^(fc') = (i^(^')'^'"'r(2 - d/2) . (4.5) 
For the rest of this section, we shall mainly concentrate ourselves on the d = 2 case. We are thus interested in 

1 dx 

Fd=2{k^) = :r- / -r. ^7^^ ?\ ? ■ (4-6) 

47r Jq x[1 — x)k'^ + x(m^ — mjj + 
Setting t = k'^ and dropping the irrelevant prefactor, we may focus on 

dx 

F{t) = / — -p, ^ 2^ 9 ■ (4-7) 

Jq x[1 — x)k^ + x[mf — TO2) + ™2 

Subsequently, 

dl'+^F{t) = /' , ^ = 1)! /' — ^ , , 

* ^' * J„ x{l-x)t + a ^ ' ^ ' J„ a;(l - x) (t + a)"+2 



(4.8) 



where we temporarily set a = a{x) — ™2)+™2 ^ j)oing so, we have 



^(^) = / li- ^ — 77—^ ) ■ (4-9) 



I ^ — 7 1™ — i^^r 7 — 7 

/q x(1 — n->-+oo n\ \[t + a) 
We already computed the derivative appearing in the r.h.s. of (|4.9|l, reusing the result p.8[) yields 



n->-+oo (n!)2 _/q x[1 — x) (t + a)^"+2 n->-+oo (n!)"' Jq x{1 — x) a 

(4.10) 

We notice the great power of this formulation, as the spectral density can now be obtained from a steepest descent 
evaluation of the Feynman parameter integral ( |4.10| i. It is useful to perform the substitution x — to rewrite 
( pTTOl i as ^ 

«^+oo (n!)2 {l + y){ml+yml) 

with 

a = a{y)^ -{jnl+yml). (4.12) 

y 
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We are now ready to search for the maxima of h{y, t) — h{a{y),t) for y G [0, oo]. Solving dyh{y, t) — gives 



mi 



t — — ml + i/— 4m^m2 + (— t + mf + m^Y 



2 



mi t — 771? — mo — \J —Ari{\mk + (— i + m? + mn)^ 

ya = , ?/4 = ^-^ • (4.13) 

m2 zmj 

We can discriminate between a few cases. 

1. < < < (mi - miY 

In this case, the root appearing in ( |4.13[ ) exists. However, only y\ > 0. We have 

9>\ - 2m3(t^ (mi+m,)^) ^ < ^ , 

« 'v-yi mi{mi + m2) (t + (mi + 1712) ) 



We can then look at the behaviour of the integral for n +oo, 

_^{n+l)h{a,t) ^ p{n+l)h{yi,t) ^ ^4 -^^-j 



_ g(«+l)/i(a,t) ^ ^ 

(1 + y){ml + yml) + T 



Using 



we find 



M.i,^)=ln JP-^^, (4.16) 
+ (mi + m2) ) 



Jn+l)h(a,t) „ 1 / t(mi+m2) ^ 



g^n+lJn^a,tJ ^ ^ . (4.17) 



(l + y)(m2 + ym2) ^/n + l \{t + [mi + m2f Y 

t{mi+m2)^ 1 
(t+(mi+m2)2)2 4 



Using (1TT61), we conclude that p{t) = as < ul^Z'^ZlW < 



2. (mi - m2)^ <t < (mi + m2)^ 
In this interval, the roots appearing in ( |4. 13| l are not real, so again only yi is relevant. The rest of the 
reasoning remains valid, so we conclude again that p{t) = 0. 

3. t > (mi + m2)^ 

In this case, the roots are again real-valued. Looking at ( |4.18| l, we notice that yi now corresponds to a 
minimum, so we can disregard it. Only j/2 and j/4 are relevant, as both are positive and correspond to local 
maxima. We shall not write down some intermediate results, as they are quite lengthy. We set 

dlh\ =-/3|<0, = dlh\ =-/32<o. (4.18) 
It then suffices to notice that 

hiy2,t)^hiy3,t) = -2\n2, (4.19) 
which is exactly the value we need to find a nonvanishing finite contribution to the integral. In particular, 

1 1 \^ 

^^^^ " (l + 2/2)(m2+y2m2) 2p2 ^ [l + yi){ml + y.m^) Wi ' ^^'^^^ 



which simplifies to 



after a bit of algebra. 
At the end of the day, we find that 



^ —Am\m\ + [t — m\ — m\Y 



Fd^2{k') = / dt-^ , with Pit) = , . , . . (4.22) 

It is easily verified that p{t) > for t > (mi + m2)^. These results are in full agreement with those obtained from 
the Cutkosky rules, see also lUl, given that we transform from Minkowski to Euclidean space. 
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5 A comment about divergent spectral integrals 



Setting d — Ain ( |4.5| l, we clearly stumble upon a divergence cx r(0). In order to deal with finite quantities, one 
usually turns to a suitably subtracted spectral representation. In particular, if we formally set 



oo 



G(x)= / dt^(=oo) (5.1) 
Jo t + x 

for a divergenj^ quantity G{x), it is clear that a suitable number of derivatives w.r.t. x will render us with a finite 
result, 

leading to a finite subtracted spectral representation after N integrations from to x, 

O'M - - ^ . . - ^I^(O) ^ (-1)"-" f "ij^^^ < - . (53, 

In order to know the spectral representation corresponding to G{x), it looks preferential to study the finite function 
H{x), (|5.2|i, and bring it into a spectral representation of the form 



oo 



from which we can also read off the desired spectral density p{t) of the original function G{x). Using the tech- 
niques of Section 2, it is not difficult to show thaj^ 

p{t) = lim [f^+^dr'^^'Hit)] , (5.5) 

a formula which allows to compute the spectral density p{t) from knowledge of the quantity H{x), and this without 
encountering any infinities at any time. 



In practice, glancing again at ( |4.5| l, we would have 



and upon comparing with ( |4.6[ ), we notice that these are almost identical. Only the "prefactor function" of the 
will be different. In particular, this means that the saddle point equation etc for d — A remains the same as for 
d = 2, just as almost all the rest of the derivation. Only during the last step in computing the spectral density, we 
shall notice a difference, as only there the "prefactor function" plays a role. 



The case d = 3 would need to be worked out from the beginning, but similar tricks as for d 2 or d = 4 can be 
applied. This can be easily appreciated a 
form that will appear in the d = 3 case. 



applied. This can be easily appreciated as a systematics in the derivatives of -yf^ exists, the latter being the basic 



6 A second application: two complex conjugate masses 

We now come to the most interesting application of this paper. We reconsider the integral ( |2.7[ ), rewritten in the 
form ( |4.6| l, but now we immediately assume that the occurring mass scales are complex conjugate, and given by 
( |2.8| l. It can be checked that the "trick-inspired" approaches of |36 | to construct the spectral density, do not longer 
work out in the case that b ^ 0. We wish to confirm that the results ( |2.16| l and ( |2.19| l obtained via a blind trust in 
the Cutkosky approach, are actually correct. 



3. We assume here that the divergency is due to the large t-behaviour of p{t), something which is usually the case. 

4. We may always take n > N. 
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We start our analysis from the first line of ( |4.10| i. We have 

{2n + iy. dx a''t'^+^ 



ra-f+oo Jq x(1 ~ X) (t + a) 



2ri+2 



now with 



, , xim^ — mo) + mo a 2x — 1 

a = a{x) = ^ — = — r + ib- 



(6.1) 



(6.2) 



x{l — x) x{l — x) x{l — x) 

In principle, one might try to do a steepest descent evaluation of the integral \6.\\ , keeping in mind that this time 
one would need to deform the integration contour to follow a line of steepest descent through the equivalent of the 
saddle points ( |4.13| l in the complex a; -plane. This looks as a very complicated exercise, instead we shall follow a 
somewhat different reasoning. Let us first notice that \6.l\ defines a real function, as 



dx 



nj-n+l 



x{l -x){t + a)2«+2 



dx 1 



x{l — x) a \{t + a) 



at 



dx 



x{l-x) 



at 



n+l 



n+l 



dx 1 



at 



a \ {t + aY 



1/2 x{l — x) a \ (t + a) 
1/2 dx 



x{l-x) 



n+l 



a*t 



a* \{t + a*Y 



n+l 



via the substitution a: — 1 — x in the 2nd integral. We can thus say 

p{t) = 2Re[i?(t)], 

and it suffices to study 

(2n + l)! /"i/^ dx 1 



R{t) = lim 



+00 (n\) 



02 



.t(1 — x) a \{t + a) 



at 



n+l 



(6.3) 



(6.4) 



It would be beneficial if we could use a as (complex) integration variable, whereby we integrate along the contou:[^ 
7, which has the parameter representation 



7 : a; e [0, 1/2] a{x) 



ib- 



2x - 1 



x{l — x) x{l — x) 



(6.5) 



This 7 will start from a point at complex infinity in the right lower half plane and will end on the real axis at 
a = 4a. Inverting gives 



-2ib + a± - iaa + 

2a 



(6.6) 



A careful examination learns that the minus sign is the appropriate choice, given that a{x) lives in the lower half 
a-plane for x G [0, 1/2]. We can subsequently reexpress (|6.4|i as 



(2n+l)! 
i? t = lim ^ , ,J 



-de 



1 



y/-462 + a(-4a + a) a \{t + a) 



at 



n+l 



(6.7) 



after a little algebra. For the further analysis, it is easier to switch the orientation of the contour 7, so than we can 
write 



R{t) = lim 



n—>-\-oo 



(2n + 1)! 
(n!)2 



da 



1 



,(n+l) In h(a,t) 



^-462 + a(-4a + a) a 



(6.8) 



with r as shown in Figure 1, and h{a, t) given by p.lO| i. 



5. We do not necessarily refer to a closed contour here! 
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Figure 1 



The presence of the square root in ( |6.8| l compUcates the situation a bit. Setting a = a^; + iay learns there is a 
"crucifix-cut" given by 

= 2a , ay € R 

ax G [2(a - Va^ + b^), 2{a + + b^)] , = . ^ 



Our goal is now to compute ( |6.8| l by means of a steepest descent approach. We can use previously gained knowl- 
edge, from which we learn that h{a,t) has a saddle point at a = t. In the vicinity of the saddle point, we may 
write ^ 

ln/t(Q;,t) = -21n2- +.... (6.10) 

We can now consider a few cases. 

1. t>2{a + Va'^ + b^) 

Since everything is analytic in the considered region of the complex plan^ we can deform our contour F 
into r', which shares its begin and end point with F, such that F' passes through the saddle point. We have to 
choose the orientation of the contour F' in the vicinity of the saddle point in such a way that the imaginary 
part of In h remains constant. From ( |6.10[ ), it is easily seen that this is the case when we take a real for 
a while until we have passed t, and then we can let F' bend over to let it flow to the correct end point at 
infinitjj^ Doing so, we find 

lim f '^^ lp(«+l)ln/.(a,f) 



r ^-4&2 -(- a(-4a + a) a 

lim ^ 1 r°°da.e-("+i) K - 

^-Ab^ +t{-Aa + t) t At 

lim ^ (6.11) 

n^oo 4«+i yrr^^ ^-462 + t{-Aa + t) 



where we made use of p.l4[ ). Combining this result with p.l6| l gives 

1 

V-452 + t(-4a + t) 



R{t) = (6.12) 



Hence, 



p{t) = , (6.13) 
V-462+i(-4a + t) 



6. The presence of the In is of no concern for the analyticity, since n S N. 

7. As the function vanishes at infinity, we can also move the end point at infinity, if desired. 
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4a < t < 2(a + Va^ + &2) 

In this region, we can still deform the contour F into an appropriate F' in a completely similar way, as we do 
not have to cross any cut in order to do so. Since we are below the cut, we find after a completely analogous 
reasoning as in the first case 



^462 + i(4a - t) 



so that now 



Pit) 



0. 



(6.14) 



(6.15) 



2a <t <4a 

At first sight, we can repeat the analysis of subcase 2. This reasoning is however flawed. A basic ingredient 
of the steepest descent approach is that Re(ln h{a, t)) reaches its maximum along the (deformed) contour at 
the value of the saddle point, in this case a — t.lf we have to first stretch our contour F into the left direction 
to pick up the saddle point, we are violating this assumption, as it can be checked. The grey region in Figure 
2 displays the region in the complex a-plane we cannot cross with our contour. 





/ 2 
/ 1 
1 . . . 1 . . . 


W 


-■6 


-4 -J ■ 

\ 

V —2 





Figure 2 



The boundaries of these curves are defined by 



tty = ± J 7<2 - 4:^yp{3t-2a.J;) - 2ta^ 



yI, ay=±\l 7t2 - 2ta^ -al + 'i^M'^ ~ 2t^a^ . 

(6.16) 

Mt > 4a, the begin point of F lies within the white blob, and we can (only) follow the a^j-axis, pick up the 
saddle point at ax = t, leave this region and go to infinity. If however we start outside this region, there is 
no way that we can go through the saddle point and go back to infinity without self-crossing the contour. 
Said otherwise, the steepest descent approach is inapplicable in this case. 

< t < 2a 

In this region, we seem to be in even more serious trouble, as we can even no longer deform our contour in 
a useful way, due to the vertical branch cut at a = 2a. 

< t < 4a 

Let us now present a combined study of the two previous cases, based on a different approach. Using the 
triangle inequality for integrals, we have 



\R{t)\ 



< 



lim 



lim 

n— f +00 



lim 

n— >+oo 



(2n + l)! 


/•1/2 


(n!)2 


Jo 


(2n + l)! 


.1/2 


(niy j 





(2n + 1)! 


.1/2 


{niy j 






dx 



x{l-x) (t + a)2»+2 



dx 
x{l — x) 

dx 



{t + a)2"+2 

1 



tv/a2 + 62(i_2a;)2a;(l - x) 
;(l-a;) H I (a - - .t))2 + &2(i _ 2.t)2 



n+l 



. (6.17) 
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We consequently observe that the (positive) function 

+ fe2(i _ 2x)^x{l - x) 



(a-tT(l-a;))2 + 62(l-2a;)2 



has maxima at the solutions of 



namely, 
with 



{x,t) = (0,0), 





dj 


dx 


dt 


{x,t) 


= (1,0), 



j(l,0) = 0, j(l/2,4a) = 



1 



j(0,0)-0. 

On the interval (x, t) e [0, 1/2] x [0, 2a], this implies that, for a certain C, 

0<j(a;,t)<C< ^, 
and using this we can majorate ( |6.17| i further as 

\R{t)\ < lim \ / C"+' / -r, -^ri 

^ ^ 

number independent of n 



Once more using p.l6| >, we then simply find 

\m\ < 0, 

as ( < i. A fortiori, R{t) — 0, and thus also p{t) = 0. 
Combining all information gathered so far, we have shown that 



r+oo 
J2(a+\/a 



dt 



Pit) 



with p{t) 



2(a+\/a2+62) t + ' 

whereby we can observe that p{t) > for t > 2{a + ^/a^~+lp). 



27r^-462 + t(-4a + <) 



(6.18) 

(6.19) 
(6.20) 
(6.21) 

(6.22) 
(6.23) 



(6.24) 



(6.25) 



As a check on the final result ( |6.25| l, we notice that its & ^ limit coincides with the -> ni^ —t' a limit of the 
earlier obtained result ( |4.22| i. 



For the general case, we could of course explicitly compute the spectral integral of ( |6.25[ ), as well as the original 
Feynman parameter integral ( |4.6| l, and verify if both results are the same. We find it however more instructive for 
the reader to display in Figure 3 both results for the explicit example ml = 1 + 2i, = 1 — 2i, making their 
equivalence clearly visible. 



(a) 



Figures: Plots of Fi{k'^ 



dx 



Q x(l — x)t + 2ibx + a — ib 



andF2(fc^ 



(b) 



rft-^^fora= 1,6 = 2. 



13 



Similar conclusions can be reached for d = 3 or d = 4. In all cases, we thus find perfect agreement between 
the results obtained by application of the Cutkosky cut rules, and those obtained by using the Widder Stieltjes 
formalism. 



7 Concluding remarks 

We have demonstrated the usefulness of the Stieltjes integral transform to compute the Kallen-Lehmann spectral 
density of (Euclidean) two-point functions, even in cases where the validity of the usually employed Cutkosky 
rules is not clear The main recipes are given in Sections 4 and 5, in particular in expression p.2| i which gives the 
spectral density p{t) in terms of the two-point function F{k'^), when the latter is known for Euclidean momenta 
> 0. 

We have paid particular attention to the case of Gribov-hke propagators, which entails the presence of propagators 
with 2 complex conjugate masses. It turns out that, at least at one loop, this case can we worked out analytically 
to the end using a steepest descent approach. Our main conclusion is that the results, obtained by the cut rules for 
real masses and in Minkowski space, may be continued to Euclidean space, while also the two real masses can 
formally be replaced by two complex conjugate values. This kind of results already found use in recent works like 
and have now been proven to be correct. 



However, the inversion formula ( |3.2| i is also valid in cases where the cut rules are unavailable. We recall that many 
works have been devoted to the study of the quantum equations of motion (Schwinger-Dyson formalism), giving 
also estimates for the gluon propagator, amongst other quantities 191 ITOl [TSl [T6l [TTl |2T|| . In the deep infrared and 
ultraviolet, an analytical approach is possible. In the intermediate momentum regime, however, numerical solutions 
are in order due to the complexity of the system. The numerical outcome can then fitted by well-guessed functional 
forms. Let us mention two major examples 



1 TTl 



4 15 

f{k\ml) = piml + p2 ^ , Pi = -x, P2 = x, TOo = 612MeV, A = 645MeV , (7.1) 

which represents the so-called "massive" ifTSi (also known as "decoupling" ifTTl ) solution, or 

13 

r a„ uz / ^ \-z \ -\ — 

D{k^) = Zijii^ 



a(0) 4tt k'^ f 1 A 



l + fc2/A2 l3ok^ + A^\\nk^/A^ P _ ^2 
.91 

'n 



/3o = a(0) = ^, A = 710MeV, (7.2) 



with 



Ziik') = wjjj^^^^^^^, orZnie) = wn^-j^^^j , 0.595 , (7.3) 

for the so-called "scaling" solution B3l . 

It is evidently out of the question to apply the standard cut rules with the foregoing propagators, since they are not 
even of the form p2^„^2 , so the replacement rule p.3| l loses its meaning. The theory of the (Stieltjes) inversion 
could however be applied, albeit perhaps rather in a numerical fashion, employing stable approximations for the 



infinite numbers of derivatives appearing in formula p.2| i. Let us also mention here that in many cases, the concrete 
application of the Cutkosky rules become intractable beyond one loop, due to the highly complicated phase space 
integrals. 

In the particular case of bound states, like a glueball or meson, which can crucially depend on a viable input for the 
gluon propagator, one might imagine to use propagators of the type ( |7.1| i or ( |7.2| i to compute, in one approximation 
or another, suitable two-point functions, to obtain consequently an estimate for the spectral density, and then 
eventually use the latter as input for whatever method one likes to employ to find estimates for e.g. bound state 
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masses. One might think about employing Laplacian sum rules as in ||6l or techniques based on Pade approximants 
as in B0ll45]| . to name only a few approaches to the bound state problem whereby the spectral density enters. 

We end by mentioning that ( |7.1| i or ( |7.2| i are quite different in the deep infrared as ( |7.1| i tends to a strictly positive 
constant for 0, while ( |7.2| l vanishes in the same limit. The former scenario is recovered on the lattice 

lfT2l [131 [T4l . One could nevertheless wonder if both solutions could reproduce more or less the same physics, as 
preliminary investigated in |46| in a specific example and approximation. We hope to come back to this issue in 
the future, using the tools developed in this paper The content of Section 5 should be of particular interest, as also 
the propagators ( |7.1| i or ( |7.2| l will lead to divergent spectral integrals. 
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